DNA sequence recognition by hybridization to short oligomers: experimental verification of the method on the E. coli genome.

A newly developed method for sequence recognition by hybridization to short oligomers is verified for the first time in genome-scale experiments. The experiments involved hybridization of 15,328 randomly selected 2-kb genomic clones of Escherichia coli with 997 short oligomer probes to detect complementary oligomers within the clones. Lists of oligomers detected within individual clones were compiled into a database. The database was then searched using known E. coli sequences as queries. The goal was to recognize the clones that are identical or similar to the query sequences. A total of 76 putative recognitions were tested in two separate but complementary recognition experiments. The results indicate high specificity of recognition. Current and prospective applications of this novel method are discussed.

Discovering distinct genes represented in 29,570 clones from infant brain cDNA libraries by applying sequencing by hybridization methodology.

To discover all distinct human genes and to determine their patterns of expression across different cell types, developmental stages, and physiological conditions, a procedure is needed for fast, mutual comparison of hundreds of thousands (and perhaps millions) of clones from cDNA libraries, as well as their comparison against data bases of sequenced DNA. In a pilot study, 29,570 clones in duplicate from both original and normalized, directional, infant brain cDNA libraries were hybridized with 107-215 heptamer oligonucleotide probes to obtain oligonucleotide sequence signatures (OSSs). The OSSs were compared and clustered based on mutual similarity into 16,741 clusters, each corresponding to a distinct cDNA. A number of distinct cDNAs were successfully recognized by matching their 107-probe OSSs against GenBank entries, indicating the possibility of sequence recognition with only a few hundred randomly chosen oligomers.

Genome-scale DNA sequence recognition by hybridization to short oligomers.

Recently developed hybridization technology (Drmanac et al. 1994) enables economical large-scale detection of short oligomers within DNA fragments. The newly developed recognition method (Milosavljevic 1995b) enables comparison of lists of oligomers detected within DNA fragments against known DNA sequences. We here describe an experiment involving a set of 4,513 distinct genomic E.coli clones of average length 2kb, each hybridized with 636 randomly selected short oligomer probes. High hybridization signal with a particular probe was used as an indication of the presence of a complementary oligomer in the particular clone. For each clone, a list of oligomers with highest hybridization signals was compiled. The database consisting of 4,513 oligomer lists was then searched using known E.coli sequences as queries in an attempt to identify the clones that match the query sequence. Out of a total of 11 clones that were recognized at highest significance level by our method, 8 were single-pass sequenced from both ends. The single-pass sequenced ends were then compared against the query sequences. The sequence comparisons confirmed 7 out of the total of 8 examined recognitions. This experiment represents the first successful example of genome-scale sequence recognition based on hybridization data.

Clone clustering by hybridization.

DNA sequencing by hybridization (SBH) Format 1 technique is based on experiments in which thousands of short oligomers are consecutively hybridized with dense arrays of clones. In this paper we present the description of a method for obtaining hybridization signatures for individual clones that guarantees reproducibility despite a wide range of variations in experimental circumstances, a sensitive method for signature comparison at prespecified significance levels, and a clustering algorithm that correctly identifies clusters of significantly similar signatures. The methods and the algorithm have been verified experimentally on a control set of 422 signatures that originate from 9 distinct clones of known sequence. Experiments indicate that only 30 to 50 oligomer probes suffice for correct clustering. This information about the identity of clones can be used to guide both genomic and cDNA sequencing by SBH or by standard gel-based methods.

Sequencing by hybridization: towards an automated sequencing of one million M13 clones arrayed on membranes.

An immediately applicable variant of the sequencing by hybridization (SBH) method is under development with the capacity to determine up to 100 million base pairs per year. The proposed method comprises six steps: (i) arraying genomic or cDNA M13 clones in 864-well plates (wells of 2 mm); (ii) preparation of DNA samples for spotting by growth of the M13 clones or by polymerase chain reaction (PCR) of the inserts using standard 96-well plates, or plates having as many as 864 correspondingly smaller wells; (iii) robotic spotting of 13,824 samples on an 8 x 12 cm nylon membrane, or correspondingly more, on up to 6 times larger filters, by offset printing with a 96 or 864 0.4 mm pin device; (iv) hybridization of dotted samples with 200-2000 32P-labeled probes comprising 16-256 10-mers having a common 8-mer, 7-mer, or 6-mer in the middle (20 probes per day, each hybridized with 250,000 dots); (v) scoring hybridization signals of 5 million sample-probe pairs per day using storage phosphor plates; and (vi) computing clone order and partial-to-complete DNA sequences using various heuristic algorithms. Genome sequencing based on a combination of this method and gel sequencing techniques may be significantly more economical than gel methods alone.

DNA sequencing by hybridization: 100 bases read by a non-gel-based method.

Determination of the sequences of human and other complex genomes requires much faster and less expensive sequencing processes than the methods in use today. Sequencing by hybridization is potentially such a process. In this paper we present hybridization data sufficient to accurately read a known sequence of 100 base pairs. In independent reactions, octamer and nonamer oligonucleotides derived from the sequence hybridized more strongly to this DNA than to controls. The 93 consecutive overlapping probes were derived from a 100-base-pair segment of test DNA and additional probes were generated by incorporation of a noncomplementary base at one of the ends of 12 of the basic probes. These 12 additional probes also had a full-match target in one of the control DNAs. The test and one of five control DNAs spotted on nylon filters were hybridized with 83 octamers and 22 nonamers under low-temperature conditions. A stronger signal in DNA containing a full-match target compared to DNA with only mismatched targets was obtained with all 105 probes. In 3 cases (2.9%), the difference of signals was not significant (less than 2-fold) due to inefficient hybridization and the consequently higher influence of background. The hybridization pattern obtained enabled us to resequence the 100 base pairs by applying an algorithm that tolerates an error rate much higher than was observed in the experiment. With this result, the technological components of large-scale DNA sequencing using the sequencing by hybridization method are in place.

An algorithm for the DNA sequence generation from k-tuple word contents of the minimal number of random fragments.

An algorithm is described for generation of the long sequence written in a four letter alphabet from the constituent k-tuple words in the minimal number of separate, randomly defined fragments of the starting sequence. It is primarily intended for use in sequencing by hybridization (SBH) process- a potential method for sequencing human genome DNA (Drmanac et al., Genomics 4, pp. 114-128, 1989). The algorithm is based on the formerly defined rules and informative entities of the linear sequence. The algorithm requires neither knowledge on the number of appearances of a given k-tuple in sequence fragments, nor the information on which k-tuple words are on the ends of a fragment. It operates with the mixed content of k-tuples of the various lengths. The concept of the algorithm enables operations with the k-tuple sets containing false positive and false negative k-tuples. The content of the false k-tuples primarily affects the completeness of the generated sequence, and its correctness in the specific cases only. The algorithm can be used for the optimization of SBH parameters in the simulation experiments, as well as for the sequence generation in the real SBH experiments on the genomic DNA.

Deletion analysis of Duchenne muscular dystrophy using cDNA probes and multiplex PCR.

The authors studied Duchenne muscular dystrophy (DMD) in 14 Yugoslav families with 16 diseased men. All the sixteen patients showed typical phenotypic features of DMD. The linkage analysis was done with five probes 754, C7, PERT87.15, pERT87.30, and pERTBir. A deletion with pERT87.15 (DX 164) in one family was found. Deletion analysis was mainly performed using cDNA probes and multiplex PCR analysis. Preferential deletions were in the hot spot deletion region covered by cDNA probe 8. Among 14 examined patients with cDNA hybridization deletions were found in 5 cases (5/14). For quick and prenatal analysis the authors used multiplex PCR. With this methods deletions were found in 4 cases (4/15). The percentage of the observed deletions cases was 50% (8/16) in Yugoslav population, which shows good correlation with other European and American populations studied. Prenatal diagnosis was done for two consultants with highly ambiguous carrier risk, and no deletions were found.

Reliable hybridization of oligonucleotides as short as six nucleotides.

Although there are many new applications for hybridizing short, synthetic oligonucleotide probes to DNA, such applications have not included determining unknown sequences of DNA. The lack of clear discrimination in hybridization of oligo probes shorter than 11 nucleotides and the lack of a theoretical understanding of factors influencing hybridization of short oligos have hampered the development of their use. We have found conditions for reliable hybridization of oligonucleotides as short as seven nucleotides to cloned DNA or to oligonucleotides attached to filters. Low-temperature hybridization and washing conditions, in contrast to the high stringency conditions currently used in hybridization experiments, have the potential for allowing the simple use of all oligos of six nucleotides or longer in meaningful hybridizations. We also present the hybridization discrimination theory that provides the conceptual framework for understanding these results.

Origin of rat beta-globin haplotypes containing three and five genes.

We have reported in rat three adult beta-gene haplotypes containing either five or three genes. Detailed sequence analysis reveals that the leftmost gene is the major gene and that at the opposite end downstream lies the minor gene. All of the genes lying between them are minor-major hybrids indicating their origin by unequal crossing-over. In two haplotypes beta-globin genes were found with an L1(1) element inserted directly into IVS2. The described results allow the formulation of a pathway of mutational events leading from the ancient two-beta-gene rodent ancestor through a three-gene haplotype to five-gene haplotypes, one of which is postulated to have arisen in common laboratory strains since their capture in the wild.

Variant chromosomal arrangement of adult beta-globin genes in rat.

The genomic organization of three haplotypes of beta-globin genes was determined to resolve the question of the number of those genes in rat. Haplotype a, found in inbred strain DA, has three genes or pseudogenes, while haplotypes b, found in AO, Y5 and Wistar strains, and c, found in Wistar strain, have five genes or pseudogenes each. In haplotypes b and c, the first gene is of beta major type and the remaining four are of beta minor type. Partial sequencing of six out of 13 genes shows that duplications of beta minor genes are causing polymorphism in a number of genes. Also, in haplotype b two beta minor genes have a 6.5-kb intron 2, while in haplotype c only one beta minor gene contains such a large intron 2. The three structurally different haplotypes described are not interconvertible by single recombination events. The results indicate that the rat has the highest number of adult beta-globin genes found in mammals so far.

Sequencing of megabase plus DNA by hybridization: theory of the method.

A mismatch-free hybridization of oligonucleotides containing from 11 to 20 monomers to unknown DNA represents, in essence, a sequencing of a complementary target. Realizing this, we have used probability calculations and, in part, computer simulations to estimate the types and numbers of oligonucleotides that would have to be synthesized in order to sequence a megabase plus segment of DNA. We estimate that 95,000 specific mixes of 11-mers, mainly of the 5'(A,T,C,G)(A,T,C,G)N8(A,T,C,G)3' type, hybridized consecutively to dot blots of cloned genomic DNA fragments would provide primary data for the sequence assembly. An optimal mixture of representative libraries in M13 vector, having inserts of (i) 7 kb, (ii) 0.5 kb genomic fragments randomly ligated in up to 10-kb inserts, and (iii) tandem "jumping" fragments 100 kb apart in the genome, will be needed. To sequence each million base pairs of DNA, one would need hybridization data from about 2100 separate hybridization sample dots. Inevitable gaps and uncertainties in alignment of sequenced fragments arising from nonrandom or repetitive sequence organization of complex genomes and difficulties in cloning "poisonous" sequences in Escherichia coli, inherent to large sequencing by any method, have been considered and minimized by choice of libraries and number of subclones used for hybridization. Because it is based on simpler biochemical procedures, our method is inherently easier to automate than existing sequencing methods. The sequence can be derived from simple primary data only by extensive computing. Phased experimental tests and computer simulations increasing in complexity are needed before accurate estimates can be made in terms of cost and speed of sequencing by the new approach. Nevertheless, sequencing by hybridization should show advantages over existing methods because of the inherent redundancy and parallelism in its data gathering.

Limited polymorphism of both classes of MHC genes in four different species of the Balkan mole rat.

We analyzed the restriction fragment length polymorphism of class I and class II MHC genes in DNA from 20 individuals belonging to the four different species of the complex of species of Balkan mole rats Spalax leucodon captured at four different localities in Yugoslavia. All populations were tested with four restriction enzymes and one conserved mouse probe for each of the two classes of MHC genes. The probes employed detect either limited polymorphism of class I genes or lack of polymorphic bands containing class II genes. Of the two other subterranean rodents that have been studied, four karyotype forms of the Israeli mole rat show polymorphism in both classes of MHC genes similar to the one found in all other mammals (Nizetic et al. 1985), and the Syrian hamster shows limited polymorphism of class I genes and high polymorphism of class II genes (McGuire et al. 1985). Balkan mole rats belong to a new group in this respect, different from all mammals studied so far, since they apparently show limited polymorphism of both classes of MHC genes.

Identification of a testis-specific gene from the mouse t-complex next to a CpG-rich island.

We have used an approach based on the observation of CpG-rich regions near the 5' end of many genes to screen a panel of cosmids derived from the t-complex and tested candidate sequences for evidence of transcription in a number of different mouse tissues. One gene so identified is expressed specifically in testicular germ cells and maps to a subregion of the t-complex also containing loci involved in transmission ratio distortion and male sterility. The transcript is first detected during the pachytene stage of the first meiotic division, but is expressed in highest levels in the later haploid spermatogenic stages. Sequence analysis verified the existence of a CpG-rich element on the 5' end of the gene and predicts a unique protein species with no significant homologies to those previously determined.

Characterization of two rat globin cDNA clones.

The rat globin gene system is suitable for studying a coordinated regulation of seven genes from two gene families. A rat reticulocyte cDNA globin library has been constructed and two clones analyzed in detail. pBRrg 5 contains alpha while pBRrg X contains beta type sequence. These cloned cDNAs will be useful probes of structure and function of rat globin genes. The deduced amino acid sequences extend the information about the variability of rat globin chains.